Negative Coupling: The Coincidence of Premating Isolating Barriers Can Reduce Reproductive Isolation.

Speciation can be mediated by a variety of reproductive barriers, and the interaction among different barriers has often been shown to enhance overall reproductive isolation, a process referred to as "coupling." Here, we analyze a population genetics model to study the establishment of linkage disequilibrium (LD) among loci involved in multiple premating barriers, an aspect that has received little theoretical attention to date. We consider a simple genetic framework underlying two distinct premating barriers, each encoded by a preference locus and its associated mating trait locus. We show that their interaction can lead to a decrease in overall reproductive isolation relative to a situation with a single barrier, a process we call "negative coupling." More specifically, in our model, negative coupling results either from sexual selection that reduces divergence at all loci, or from reduced LD that occurs because the presence of many females with "mismatched" preferences causes the mating success of recombinant males to become high. Interestingly, the latter effect may even cause LD among preference loci to become negative when recombination rates among loci are low. We conclude that coincident reproductive barriers may not necessarily reinforce each other, and that the underlying loci may not necessarily develop a positive association.

Phenotypic and genomic dissection of colour pattern variation in a reef fish radiation.

Coral reefs rank among the most diverse species assemblages on Earth. A particularly striking aspect of coral reef communities is the variety of colour patterns displayed by reef fishes. Colour pattern is known to play a central role in the ecology and evolution of reef fishes through, for example, signalling or camouflage. Nevertheless, colour pattern is a complex trait in reef fishes-actually a collection of traits-that is difficult to analyse in a quantitative and standardized way. This is the challenge that we address in this study using the hamlets (Hypoplectrus spp., Serranidae) as a model system. Our approach involves a custom underwater camera system to take orientation- and size-standardized photographs in situ, colour correction, alignment of the fish images with a combination of landmarks and Bézier curves, and principal component analysis on the colour value of each pixel of each aligned fish. This approach identifies the major colour pattern elements that contribute to phenotypic variation in the group. Furthermore, we complement the image analysis with whole-genome sequencing to run a multivariate genome-wide association study for colour pattern variation. This second layer of analysis reveals sharp association peaks along the hamlet genome for each colour pattern element and allows to characterize the phenotypic effect of the single nucleotide polymorphisms that are most strongly associated with colour pattern variation at each association peak. Our results suggest that the diversity of colour patterns displayed by the hamlets is generated by a modular genomic and phenotypic architecture.

A review of 263 years of taxonomic research on Hypoplectrus (Perciformes: Serranidae), with a redescription of Hypoplectrus affinis (Poey, 1861).

The hamlets (Hypoplectrus spp., Perciformes: Serranidae) constitute a distinctive model system for the study of a variety of ecological and evolutionary processes including the evolution and maintenance of simultaneous hermaphroditism and egg trading, sex allocation, sexual selection, social-trap, mimicry, dispersal, speciation, and adaptive radiation. Addressing such fundamental and complex processes requires a good knowledge of the taxonomy and natural history of the hamlets. Here, we review the taxonomy of the hamlets, from early ichthyological studies to the most recent species description in 2018. We report a total of 72 different binomial names for Hypoplectrus, synonymized or invalidated down to 17 unambiguously recognized species today. In addition, we redescribe Hypoplectrus affinis (Poey, 1861) as a valid species. In Bocas del Toro (Panama), this hamlet is distinct from eight sympatric congeners in terms of colour pattern, body size and behaviour. Whole-genome analysis and spawning observations indicate that it is genetically distinct from sympatric congeners and reproductively isolated through assortative pairing. Based on the colour pattern we detail in its redescription, live-fish photographs, videos, and earlier reports, H. affinis occurs in Panama, Nicaragua, Mexico, the Florida Keys, Cuba, Grand Cayman, Jamaica, the Dominican Republic, Los Roques (Venezuela), Bonaire, and Tobago. We conclude with a discussion of pending taxonomic issues in this group and the species status of the hamlets in general.

Rapid radiation in a highly diverse marine environment.

Rapid diversification is often observed when founding species invade isolated or newly formed habitats that provide ecological opportunity for adaptive radiation. However, most of the Earth's diversity arose in diverse environments where ecological opportunities appear to be more constrained. Here, we present a striking example of a rapid radiation in a highly diverse marine habitat. The hamlets, a group of reef fishes from the wider Caribbean, have radiated into a stunning diversity of color patterns but show low divergence across other ecological axes. Although the hamlet lineage is ∼26 My old, the radiation appears to have occurred within the last 10,000 generations in a burst of diversification that ranks among the fastest in fishes. As such, the hamlets provide a compelling backdrop to uncover the genomic elements associated with phenotypic diversification and an excellent opportunity to build a broader comparative framework for understanding the drivers of adaptive radiation. The analysis of 170 genomes suggests that color pattern diversity is generated by different combinations of alleles at a few large-effect loci. Such a modular genomic architecture of diversification has been documented before in Heliconius butterflies, capuchino finches, and munia finches, three other tropical radiations that took place in highly diverse and complex environments. The hamlet radiation also occurred in a context of high effective population size, which is typical of marine populations. This allows for the accumulation of new variants through mutation and the retention of ancestral genetic variation, both of which appear to be important in this radiation.

A novel metabarcoding primer pair for environmental DNA analysis of Cephalopoda (Mollusca) targeting the nuclear 18S rRNA region.

Cephalopods are pivotal components of marine food webs, but biodiversity studies are hampered by challenges to sample these agile marine molluscs. Metabarcoding of environmental DNA (eDNA) is a potentially powerful technique to study oceanic cephalopod biodiversity and distribution but has not been applied thus far. We present a novel universal primer pair for metabarcoding cephalopods from eDNA, Ceph18S (Forward: 5'-CGC GGC GCT ACA TAT TAG AC-3', Reverse: 5'-GCA CTT AAC CGA CCG TCG AC-3'). The primer pair targets the hypervariable region V2 of the nuclear 18S rRNA gene and amplifies a relatively short target sequence of approximately 200 bp in order to allow the amplification of degraded DNA. In silico tests on a reference database and empirical tests on DNA extracts from cephalopod tissue estimate that 44-66% of cephalopod species, corresponding to about 310-460 species, can be amplified and identified with this primer pair. A multi-marker approach with the novel Ceph18S and two previously published cephalopod mitochondrial 16S rRNA primer sets targeting the same region (Jarman et al. 2006 Mol. Ecol. Notes. 6, 268-271; Peters et al. 2015 Mar. Ecol. 36, 1428-1439) is estimated to amplify and identify 89% of all cephalopod species, of which an estimated 19% can only be identified by Ceph18S. All sequences obtained with Ceph18S were submitted to GenBank, resulting in new 18S rRNA sequences for 13 cephalopod taxa.

New genomic resources for three exploited Mediterranean fishes.

Extensive fishing has led to fish stock declines throughout the last decades. While clear stock identification is required for designing management schemes, stock delineation is problematic due to generally low levels of genetic structure in marine species. The development of genomic resources can help to solve this issue. Here, we present the first mitochondrial and nuclear draft genome assemblies of three economically important Mediterranean fishes, the white seabream, the striped red mullet, and the comber. The assemblies are between 613 and 785 Mbp long and contain between 27,222 and 32,375 predicted genes. They were used as references to map Restriction-site Associated DNA markers, which were developed with a single-digest approach. This approach provided between 15,710 and 21,101 Single Nucleotide Polymorphism markers per species. These genomic resources will allow uncovering subtle genetic structure, identifying stocks, assigning catches to populations and assessing connectivity. Furthermore, the annotated genomes will help to characterize adaptive divergence.

The Evolution of Egg Trading in Simultaneous Hermaphrodites.

Egg trading-whereby simultaneous hermaphrodites exchange each other's eggs for fertilization-constitutes one of the few rigorously documented and most widely cited examples of direct reciprocity among unrelated individuals. Yet how egg trading may initially invade a population of nontrading simultaneous hermaphrodites is still unresolved. Here, we address this question with an analytical model that considers mate encounter rates and costs of egg production in a population that may include traders (who provide eggs for fertilization only if their partners also have eggs to reciprocate), providers (who provide eggs regardless of whether their partners have eggs to reciprocate), and withholders (cheaters who mate only in the male role and just use their eggs to elicit egg release from traders). Our results indicate that a combination of intermediate mate encounter rates, sufficiently high costs of egg production, and a sufficiently high probability that traders detect withholders (in which case eggs are not provided) is conducive to the evolution of egg trading. Under these conditions, traders can invade-and resist invasion from-providers and withholders alike. The prediction that egg trading evolves only under these specific conditions is consistent with the rare occurrence of this mating system among simultaneous hermaphrodites.

The evolution of microendemism in a reef fish (Hypoplectrus maya).

Marine species tend to have extensive distributions, which are commonly attributed to the dispersal potential provided by planktonic larvae and the rarity of absolute barriers to dispersal in the ocean. Under this paradigm, the occurrence of marine microendemism without geographic isolation in species with planktonic larvae poses a dilemma. The recently described Maya hamlet (Hypoplectrus maya, Serranidae) is exactly such a case, being endemic to a 50-km segment of the Mesoamerican Barrier Reef System (MBRS). We use whole-genome analysis to infer the demographic history of the Maya hamlet and contrast it with the sympatric and pan-Caribbean black (H. nigricans), barred (H. puella) and butter (H. unicolor) hamlets, as well as the allopatric but phenotypically similar blue hamlet (H. gemma). We show that H. maya is indeed a distinct evolutionary lineage, with genomic signatures of inbreeding and a unique demographic history of continuous decrease in effective population size since it diverged from congeners just ~3,000 generations ago. We suggest that this case of microendemism may be driven by the combination of a narrow ecological niche and restrictive oceanographic conditions in the southern MBRS, which is consistent with the occurrence of an unusually high number of marine microendemics in this region. The restricted distribution of the Maya hamlet, its decline in both census and effective population sizes, and the degradation of its habitat place it at risk of extinction. We conclude that the evolution of marine microendemism can be a fast and dynamic process, with extinction possibly occurring before speciation is complete.

Inter-chromosomal coupling between vision and pigmentation genes during genomic divergence.

Recombination between loci underlying mate choice and ecological traits is a major evolutionary force acting against speciation with gene flow. The evolution of linkage disequilibrium between such loci is therefore a fundamental step in the origin of species. Here, we show that this process can take place in the absence of physical linkage in hamlets-a group of closely related reef fishes from the wider Caribbean that differ essentially in colour pattern and are reproductively isolated through strong visually-based assortative mating. Using full-genome analysis, we identify four narrow genomic intervals that are consistently differentiated among sympatric species in a backdrop of extremely low genomic divergence. These four intervals include genes involved in pigmentation (sox10), axial patterning (hoxc13a), photoreceptor development (casz1) and visual sensitivity (SWS and LWS opsins) that develop islands of long-distance and inter-chromosomal linkage disequilibrium as species diverge. The relatively simple genomic architecture of species differences facilitates the evolution of linkage disequilibrium in the presence of gene flow.

Long-Distance Benefits of Marine Reserves: Myth or Reality?

Long-distance (>40-km) dispersal from marine reserves is poorly documented; yet, it can provide essential benefits such as seeding fished areas or connecting marine reserves into networks. From a meta-analysis, we suggest that the spatial scale of marine connectivity is underestimated due to the limited geographic extent of sampling designs. We also found that the largest marine reserves (>1000km2) are the most isolated. These findings have important implications for the assessment of evolutionary, ecological, and socio-economic long-distance benefits of marine reserves. We conclude that existing methods to infer dispersal should consider the up-to-date genomic advances and also expand the spatial scale of sampling designs. Incorporating long-distance connectivity in conservation planning will contribute to increase the benefits of marine reserve networks.

Diet and stable isotope analyses reveal the feeding ecology of the orangeback squid Sthenoteuthis pteropus (Steenstrup 1855) (Mollusca, Ommastrephidae) in the eastern tropical Atlantic.

In the eastern tropical Atlantic, the orangeback flying squid Sthenoteuthis pteropus (Steenstrup 1855) (Cephalopoda, Ommastrephidae) is a dominant species of the epipelagic nekton community. This carnivore squid has a short lifespan and is one of the fastest-growing squids. In this study, we characterise the role of S. pteropus in the pelagic food web of the eastern tropical Atlantic by investigating its diet and the dynamics of its feeding habits throughout its ontogeny and migration. During three expeditions in the eastern tropical Atlantic in 2015, 129 specimens were caught by hand jigging. Stomach content analyses (via visual identification and DNA barcoding) were combined with stable isotope data (∂15N and ∂13C) of muscle tissue to describe diet, feeding habits and trophic ecology of S. pteropus. Additionally, stable isotope analyses of incremental samples along the squid's gladius-the chitinous spiniform structure supporting the muscles and organs-were carried out to explore possible diet shifts through ontogeny and migration. Our results show that S. pteropus preys mainly on myctophid fishes (e.g. Myctophum asperum, Myctophum nitidulum, Vinciguerria spp.), but also on other teleost species, cephalopods (e.g. Enoploteuthidae, Bolitinidae, Ommastrephidae), crustaceans and possibly on gelatinous zooplankton as well. The squid shows a highly opportunistic feeding behaviour that includes cannibalism. Our study indicates that the trophic position of S. pteropus may increase by approximately one trophic level from a mantle length of 15 cm to 47 cm. The reconstructed isotope-based feeding chronologies of the gladii revealed high intra- and inter-individual variability in the squid's trophic position and foraging area. These findings are not revealed by diet or muscle tissue stable isotope analysis. This suggests a variable and complex life history involving individual variation and migration. The role of S. pteropus in transferring energy and nutrients from lower to higher trophic levels may be underestimated and important for understanding how a changing ocean impacts food webs in the eastern Atlantic.

Population genomics of local adaptation versus speciation in coral reef fishes (Hypoplectrus spp, Serranidae).

Are the population genomic patterns underlying local adaptation and the early stages of speciation similar? Addressing this question requires a system in which (i) local adaptation and the early stages of speciation can be clearly identified and distinguished, (ii) the amount of genetic divergence driven by the two processes is similar, and (iii) comparisons can be repeated both taxonomically (for local adaptation) and geographically (for speciation). Here, we report just such a situation in the hamlets (Hypoplectrus spp), brightly colored reef fishes from the wider Caribbean. Close to 100,000 SNPs genotyped in 126 individuals from three sympatric species sampled in three repeated populations provide genome-wide levels of divergence that are comparable among allopatric populations (F st estimate = 0.0042) and sympatric species (F st estimate = 0.0038). Population genetic, clustering, and phylogenetic analyses reveal very similar patterns for local adaptation and speciation, with a large fraction of the genome undifferentiated (F st estimate ≈ 0), a very small proportion of F st outlier loci (0.05-0.07%), and remarkably few repeated outliers (1-3). Nevertheless, different loci appear to be involved in the two processes in Hypoplectrus, with only 7% of the most differentiated SNPs and outliers shared between populations and species comparisons. In particular, a tropomyosin (Tpm4) and a previously identified hox (HoxCa) locus emerge as candidate loci (repeated outliers) for local adaptation and speciation, respectively. We conclude that marine populations may be locally adapted notwithstanding shallow levels of genetic divergence, and that from a population genomic perspective, this process does not appear to differ fundamentally from the early stages of speciation.

Pairing dynamics and the origin of species.

Whether sexual selection alone can drive the evolution of assortative mating in the presence of gene flow is a long-standing question in evolutionary biology. Here, we report a role for pairing dynamics of individuals when mate choice is mutual, which is sufficient for the evolution of assortative mating by sexual selection alone in the presence of gene flow. Through behavioural observation, individual-based simulation and population genetic analysis, we evaluate the pairing dynamics of coral reef fish in the genus Hypoplectrus (Serranidae), and the role these dynamics can play for the evolution of assortative mating. When mate choice is mutual and the stability of mating pairs is critical for reproductive success, the evolution of assortative mating in the presence of gene flow is not only possible, but is also a robust evolutionary outcome.

Estimating dispersal from genetic isolation by distance in a coral reef fish (Hypoplectrus puella).

The spatial scale of dispersal in coral reef fishes eludes ecologists despite the importance of this parameter for understanding the dynamics of ecological and evolutionary processes. Genetic isolation by distance (IBD) has been used to estimate dispersal in coral reef fishes, but its application in marine systems has been limited by insufficient sampling at different spatial scales and a lack of information regarding population density. Here, we present an analysis of IBD in the barred hamlet (Hypoplectrus puella, Serranidae) at spatial scales ranging from 10 to 3200 km complemented with SCUBA surveys of population densities covering 94000 m2 of reef. We used 10 hypervariable DNA markers to genotype 854 fish from 15 locations, and our results establish that IBD in H. puella emerges at a spatial scale of 175 km and is preserved up to the regional scale (3200 km). Assuming a normal or a Laplace dispersal function, our data are consistent with mean dispersal distances in H. puella that range between 2 and 14 km. Such small mean dispersal distances is a surprising result given the three-week pelagic larval duration of H. puella and the low level of genetic structure at the Caribbean scale (Wright's fixation index, F(ST), estimate = 0.005). Our data reinforce the importance of considering population density when estimating dispersal from IBD and underscore the relevance of sampling at local scales, even when genetic structure is weak at the regional scale.

Population genetic analyses of Hypoplectrus coral reef fishes provide evidence that local processes are operating during the early stages of marine adaptive radiations.

Large-scale, spatially explicit models of adaptive radiation suggest that the spatial genetic structure within a species sampled early in the evolutionary history of an adaptive radiation might be higher than the genetic differentiation between different species formed during the same radiation over all locations. Here we test this hypothesis with a spatial population genetic analysis of Hypoplectrus coral reef fishes (Serranidae), one of the few potential cases of a recent adaptive radiation documented in the marine realm. Microsatellite analyses of Hypoplectrus puella (barred hamlet) and Hypoplectrus nigricans (black hamlet) from Belize, Panama and Barbados validate the population genetic predictions at the regional scale for H. nigricans despite the potential for high levels of gene flow between populations resulting from the 3-week planktonic larval phase of Hypoplectrus. The results are different for H. puella, which is characterized by significantly lower levels of spatial genetic structure than H. nigricans. An extensive field survey of Hypoplectrus population densities complemented by individual-based simulations shows that the higher abundance and more continuous distribution of H. puella could account for the reduced spatial genetic structure within this species. The genetic and demographic data are also consistent with the hypothesis that H. puella might represent the ancestral form of the Hypoplectrus radiation, and that H. nigricans might have evolved repeatedly from H. puella through ecological speciation. Altogether, spatial genetic analysis within and between Hypoplectrus species indicate that local processes can operate at a regional scale within recent marine adaptive radiations.

Colour pattern as a single trait driving speciation in Hypoplectrus coral reef fishes?

Theory shows that speciation in the presence of gene flow occurs only under narrow conditions. One of the most favourable scenarios for speciation with gene flow is established when a single trait is both under disruptive natural selection and used to cue assortative mating. Here, we demonstrate the potential for a single trait, colour pattern, to drive incipient speciation in the genus Hypoplectrus (Serranidae), coral reef fishes known for their striking colour polymorphism. We provide data demonstrating that sympatric Hypoplectrus colour morphs mate assortatively and are genetically distinct. Furthermore, we identify ecological conditions conducive to disruptive selection on colour pattern by presenting behavioural evidence of aggressive mimicry, whereby predatory Hypoplectrus colour morphs mimic the colour patterns of non-predatory reef fish species to increase their success approaching and attacking prey. We propose that colour-based assortative mating, combined with disruptive selection on colour pattern, is driving speciation in Hypoplectrus coral reef fishes.